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Abstract 

The production of antiprotons in pp collisions is investigated close to thresh- 
old where experimental data about the total cross section are not available. We 
analyze the latter reaction within the LUND string model for inclusive p produc- 
tion and within the framework of a one-boson exchange model for the exclusive 
reaction pp — > pppp. The application of our new results to the analysis of sub- 
threshold antiproton production in d + p and d + d collisions shows cross sections 
that are much lower than expected before. Nevertheless, the comparison of ex- 
perimental p differential cross sections from d+p and d+d is expected to provide 
valuable information about a nonnucleonic component in the deuteronwavefunc- 
tion. 
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The production of particles at energies below the free nucleon- nucleon threshold is 
one of the most promising sources of information about the properties of nuclear matter 
at high densities or about nucleon-nucleon correlations at short N — N distances 
Apart from heavy meson (rj, K + , K~ , p, u, $) production, the investigation of antipro- 
tons is of particular interest since they involve a much larger production threshold and 
can be more easily identified with magnetic spectrometers due to their large mass and 
negative charge. In fact, antiprotons have been detected at far subthreshold energies 
in both, p + A and A + A collisions 0, |3], £|, |], ||. The actual magnitude of the 
cross sections observed indicate strong in-medium modifications of the antiprotons as 
found from independent transport theoretical studies 0, || || |T(J. Although the proper 
magnitude of the p potential and its annihilation cross section in the medium is still 
a matter of debate, it is clear that the dominant production mechanism in nucleus- 
nucleus collisions proceeds via intermediate baryon resonances since the latter act as 



short-time energy reservoirs for the p production M, [LOL 11]. On the other hand the 



iVA or even AA channels play a minor role in proton-nucleus and deuteron-nucleus 
reactions since the resonances on average decay before colliding with another nucleon 
due to the much lower densities involved [12[ . 

However, the poor knowledge of the elementary production cross section pp — > p+X 
especially close to threshold leads to large ambiguities in the interpretation of the 
experimental data about subthreshold p production in p + A and A + A collisions and 
thus to sizeable uncertainties for the antiproton potential or selfenergy in the nucleus. 
It is thus mandatory to analyse the elementary production in situations where the 
kinematical conditions are more clearly defined and where in-medium potentials as well 
as intermediate pion induced production channels as well as antiproton annihilation 
can approximately be neglected. This is clearly the case for pp collisions and also quite 
well fulfilled for d + p and d + d reactions according to the analysis in fl2], |13| . 

In this paper we thus concentrate first on the inclusive process pp — > pX, which 
in comparison to the data from [14] is quite well described within the LUND string 



formation and fragmentation model fl5| at invariant energies above ^/s ~ 4.7 GeV. 
Within the same approach we then compute the cross section for the exclusive reaction 
pp — > pppp which due to kinematical reasons is equal to the inclusive cross section 
close to threshold. The latter exclusive channel, furthermore, is described by using i) 
a constant matrix element and ii) matrix elements calculated within the framework of 
the one-boson exchange model where coupling constants and formfactors are fitted to 
other related cross sections. The sensitivity of the antiproton spectra in subthreshold 
d — p, d — d collisions then is reanalyzed in particular with respect to a non-nucleonic 
component in the deuteron wave function (cf. ref. |L3], |16[). 
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We start with the inclusive p production in pp collisions and show in Fig. la) 
the available experimental data from |14| (full squares) in comparison to the results 
from the LUND string model (LSM) |L5] (open circles) as a function of the invariant 
energy above threshold. Since the description of the inclusive data within the LSM 
is quite good, we use the same approach to compute the relative fraction of events 
corresponding to the exclusive channel pp —>■ pppp which is represented in Fig. la) by 
the crosses. Whereas at higher energies the inclusive antiproton production corresponds 
to events with 3 baryons and further mesons, below about ^/s — Am « 0.4 GeV the 
dominant channel is pp — *> pppp. It should be noted that the LSM at this energy predicts 
cross sections well below the parametrization from Batko et al. (dashed line in Fig. 
lb), however, it is not clear if the LSM should provide reasonable extrapolations closer 
to threshold. We thus have to employ a microscopic model for the latter exclusive 
channel to obtain information about the cross section closer to threshold. 

For this purpose we analyze the process pp — > pppp within the framework of the 
one-boson-exchange (OBE) model according to the diagram in Fig. 2 describing the 
pp production via the off-shell production of tt°, p° and to pairs that annihilate to a pp 
pair. The general expression for this cross section (denoted by 04) can be written in 
the form: 

1 r s i f i i f i i r s i f i t 

Ua = f / ds 3 l dt 3 I dU I dso / (1) 

i28(2tt) 5 A(s, m 2 ,m 2 ) Jm* Jt~ Jt- Js~ Jt~ KJ 

I m I \ 1 2 dt 2 

\ T NN _ >PpNN {ti,t2,t 3 , s 2 ) 1/2 , — \m/o/ 7 2T 

\ 1 / 2 {s2,t 3 ,t 1 )\V 2 (s 3 ,t 3 ,m 2 ) 

where the following notations are introduced : A 1//2 (x,y, z) — (x— (y 1 ^ 2 + z 1 ^ 2 ) 2 ) 1 ^ 2 (x — 
(y 1 / 2 — z 1 / 2 ) 2 ) 1 / 2 if y 1 / 2 > and z 1 ! 2 > ; t 1 and t 3 are the transfers from the 
initial nucleon to the final one corresponding to the upper and lower verticies of the 
graph in Fig. 2, respectively, or the squares of the four- momenta of the exchanged 
mesons in the intermediate state ; t 2 is the transfer from the intermediate meson to 
the final antiproton; S2 = m^ p is the square of the effective mass of the pp pair, and 
•S3 = (Pmes. +Pn) 2 , Pmes., Pn are the four-momenta of the intermediate meson and the 
initial nucleon N; 



si = (s^ 2 -m) 2 ; 



tf = 2m 2 --L( S3 + m 2 - t 3 ){s 3 + m 2 - s 2 ) T A 1/2 (s 3 , m 2 , t 3 )X 1/2 {s 3 , m 2 , s 2 )); 
2s 3 

t f = t 1 + m 2 - -*-(s 2 + tt - t 3 )s 2 T A 1/2 (s 2 , t u t 3 )X 1/2 (s 2 , m 2 , m 2 )); 

ZS 2 
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sf = s 3 + m 2 - 7^(s 3 + m 2 - t 3 )(2m 2 - h) T A 1/2 (s 3 , m 2 , t 3 )A 1/2 (t 1; m 2 , m 2 ); 

s™ in = mai(s2,4m 2 ) 

tf = 2m 2 - —s{s + m 2 - s 3 ) q= A 1/2 (s, m 2 , m 2 )A 1/2 (s, m 2 , s 3 )). 

The matrix element Tfq^^ pp ^fq{t\,t 2 ,tj > ,s 2 ) can be calculated within the framework 
of the one-meson exchange model taking into account both pseudoscalar, scalar and 
vector mesons fL8f . Neglecting the higher order terms in t 2 /m 4 caused by the tensor 



part of the NNp-vertex, which is legitimate at not too large transfers t\ and t 3 (cf. 
[HS[]), it can be written in the form: 

| T^-^r^!, t 2 , t 3 , S 2 ) | - ^ (ti _ m 2)2 (t 3 - m 2)2 

I fii-^pp(s 2 , t 2 ) (2) 

where fa^ pp (s 2 , t 2 ) is the amplitude for the process ii — > pp and the index i stands for 
the exchanged meson (cf. Fig. 2) while Fi is the corresponding formfactor and giNN is 
the iNN coupling constant corresponding to the exchanged meson (i = tt°, p°,u). The 



coupling constants and formfactors are taken from refs. |18|, [19], £(J . 

We incoherently sum the contributions of the 7r°, p°, u;-exchange graphs (cf. Fig. 2) 
because the amplitudes for the processes ii — * pppp are not known sufficiently well and 
more reliable theoretical approaches are not available so far. Furthermore, assuming 
off-mass shell effects in the amplitudes fu^ pp (s 2 , t 2 ) to be small, they can be related to 
the differential cross sections of the reactions ii — > pp by 

I fn-+pp{ti,t2,s 2 ) | 2 = 167rA(s 2 , m 2 , m - ) — (s 2 ,t 2 ). (3) 

at 2 

The latter can be written in the form: 

= au^isMh), (4) 

ut 2 

where <Ju^ pp {s 2 ) is the cross section for the production of the pp pair in the annihilation 
of the mesons of type i (ir,p,tj, etc.). In (4) 4>(t 2 ) is a function normalized to 1 that 
determines the ^-dependence of the differential cross section dau^ pp /dt 2 . In our actual 
computation it was choosen to be of exponential form 

0(t 2 ) = C ie xp(Bt 2 ), (5) 

where the constant C\ is determined by the normalization of <j)(t 2 ) and 5^6 — 9 
GeV 2 . We note that our results for the total p cross section will be without noticable 
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sensitivity to the actual value of B in the threshold regime where the cross section is 
dominated by phase space (see below). 

In order to calculate cru->pp(s2) we address to the experimental data for the cross 



section (jp P _„ r - 7r + p o p o 21] which can be related to the cross section cr 7r - 7r +_^« p (s2) or 



cr p o p o^pp(s 2 ) using the detailed balance principle^. The s 2 -dependence of the it + it~ 
cross section now can be approximated by the following expression: 



M A (hcf 
2s 2 (s 2 - M 2 ) 2 + M 2 T 2 ) 



<W-»(« a ) = - 4m 2 /^ 2 ) 1/2 (6) 



with C 2 = 1,M = 2.07GeV, T = O.QGeV, (fie) 2 = 0.389380 (GeV 2 mb). 

For the calculation of the cross section cr p o p o^pp(s2) the functional form (6) was 
taken, too, but the factor C 2 (~ 0.25) was fitted to the 3 data points available from 



ref. pl| . The cross section cr wtc; ->pp( s 2) > furthermore, was assumed to be the same as 
for the process p°p° — > pp since there are no data available. 

Within the rather drastic approximations described above the antiproton cross sec- 
tion then only depends on the meson formfactors and meson-nucleon-nucleon coupling 
constants. In line with [18| the meson formfactors F{(t) are taken to be of monopole 
form, i.e.: 



m - (T) 



involving a cut-off parameter Aj. The actual parameters used are: A„. = 0.7 GeV/c, A p 
= 2.0 GeV/c, A w = 1.5 GeV/c, fi NN /4x = 14.7, g 2 pNN /^ = 40.8 and ^ w /4tt = 20 
^, |T9|; the values for A n ,gl NN were taken from []2D| where they were found to yield 



a good description of it — N scattering. Note, that the cut-off parameters Aj and the 
coupling constants gwN fc> r p- and to- mesons correspond to the relativistic (energy- 
independent) one-boson -exchange potential as considered in ref. |TJ| |T9" |. 



Calculating the cross section cr 4 (1) as a function of ^Js within the parameters 
specified above we obtain the solid line in Fig. la) that describes very well the cross 
section for the exclusive channel from the LSM above y/s — 4m w 0.5 GeV. The dotted 
line in Fig. la) shows the result when including only 7r° exchange which, however, 
is suppressed as compared to the vector meson exchange contributions due to the 
lower restmass and cut-off parameter. Qualitatively, a similar result was found for pp 
production from meson-meson annihilation in ref. [|22] . 



One might worry about the validity of the boson exchange model for the pp pro- 
duction close to threshold. In this respect we additionally employ a simple phase-space 



1 A similar concept has been used by Ko and Ge in ref. p2|] , where the authors study the pp 
production by irir, -q-q, pp and uu channels in a hot fireball. 
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model assuming that the p cross section is proportional to the 4-body phase-space in- 
tegral i? 4 ( v /s,m,m, m,wip) with a constant fitted to the first experimental point 
at ^fs — Am ps 1 GeV. The result of this simple approximation is displayed in Fig. 
lb) by the dotted line and practically coincides with the result from the OBE model 
up to yfs — 4m « 0.7 GeV. Thus the p cross section close to threshold appears to be 
dominated by phase space, only. Now combining the results for the antiproton cross 
section from the OBE model and LSM in their respective kinematical regimes, we can 
use these cross sections for further applications in p + A and A + A reactions. We 
note that a good fit for the inclusive p cross section is given by 

v PP ^p+x{^s) ~ Ra{Vs, to, to, m, mp) _ - p [mb] (8) 

with D = 4 * 10~ 3 and T= 3 GeV (solid line in Fig. lb). 

There is a chance that final state interactions might increase again the p yield very 
close to threshold, but experiences with rj production in pp collisions indicate that a 
respective enhancement is limited to the energy range y/s — 4m < 40 MeV [^3[ . 

Before exploring the consequences for p cross sections in subthreshold hadron- 
nucleus and nucleus-nucleus reactions employing our new 'elementary' cross section, 
we consider the reactions d + p and d + d where medium effects can approximately 
be neglected as discussed above. Recently, these reactions were studied [|13|] to extract 
some new information about the nuclear structure at short N — N distances employing 



the extrapolation from Batko et al. |L7|]. Since the latter cross section now sizeably 

overestimates our new results by more than an order of magnitude close to threshold 

(cf. Fig. lb) the sensitivity of the Lorentz invariant p cross section due a non-nucleonic 

component in the deuteron wavefunction has to be reexamined. 

Our model for antiproton production in d + p and d + d is described in detail in 

ref. [^] and doesn't have to be repeated here. The only modification introduced is 

to replace the parametrized form of the elementary p cross section by our new results 

(8). The deuteron wavefunction (d.w.f.) employed is that obtained from the Paris 

potential |24| transformed to a relativistic version that only depends on the light cone 

variable x and the transverse momentum k t (cf. [Jl6)). The results of our calculation 

(with the Paris d.w.f.) for the Lorentz invariant cross section Epd 3 'a / 'dp^(\fs) at 0° in 

the laboratory system for d + p at a deuteron momentum of 10 GeV/c and d + d at 7 

GeV/c are displayed in Fig. 3 within the parametrization from Batko et al. |I7| (solid 

lines) and our new cross section (dashed lines), respectively^. Whereas the reduction 

of the cross section in the (d + p) case with the new cross section is already about a 

2 In for the calculation of the p spectra in the dp — > pX reaction a factor Att was missing in the 
normalization of the d.w.f. 
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factor of 7, the (d + d ) cross section decreases up to a factor of 36. The maximum in 
the differential cross section for (d + d) of about 7 pb c 3 /GeV 2 at P d = 7 GeV/c now 
will be hard to measure experimentally. 



We follow ref. [13] and additionally consider the possibility that the deuteron has 



a 3% admixture of a non-nucleonic component as described by Eqs. (5) - (8) in [13| in 
line with refs. |TR ^| . The result for the Lorentz invariant p cross section in this case 
is also shown in Fig. 3 using the extrapolation from Batko et al. |17j] (dot-dashed lines) 
and our new 'elementary' cross section (dotted lines). Here we find a very pronounced 
enhancement of the p yield when including the non-nucleonic component; the relative 
enhancement is even larger for the new cross section than for the parametization used 
previously. Thus by measuring antiproton production in pp collisions and comparing 
relative to the d + p reaction the existence of a non-nucleonic component or its relative 
strength should be clearly measurable. 

In ref. [[13|] it was, furthermore, suggested that the ratio of the antiproton cross 
section from d+d to d+p reactions might provide some information on the non-nucleonic 
component of the d.w.f. itself because ratios of cross sections are less sensitive to the 
actual magnitude of the elementary cross section. In fact, within the parametrization 
from ref. [13[ a relative sensitivity up to a factor of 1.5 - 2 has been found (cf. Fig. 
3 of |13|)- Our reanalysis of this suggestion with the new 'elementary' cross section is 
shown in Fig. 4 for a deuteron momentum of 9.5 GeV/c. Here the solid line reflects a 
calculation including a 3% admixture of the non-nucleonic component while the dashed 
line is obtained with the Paris d.w.f., only. Contrary to ref. [|13[ we find that the relative 
ratio dd/dp changes only slightly with the p momentum such that the ratio itself does 
no longer qualify for determining the deuteron structure. 

In summary, we have reexamined the production of antiprotons in pp, dp and dd 
reactions close to threshold energies. Our results are based on a combined analysis 
within the LUND string model |0| , an effective OBE model for the exclusive channel as 
well as on 4-body phase space and clearly indicate that the estimates for p production at 
subtreshold energies in p + A and A + A collisions within the extrapolation of Batko et 



al. [|I7| are severely overestimated. In d+p reactions at 10 GeV/c the relative reduction 
is about a factor of 7 whereas in d + d collisions at 7 GeV/c we find a reduction by a 
factor of about 36 as compared to previous estimates. These reduced cross sections, on 
the other hand, will require much more attractive antiproton selfenergies in the nuclear 
medium than estimated before in p + A and A + A reactions. 

The comparison of pp, dp and dd collisions will, however, still provide valuable infor- 
mation about a non-nucleonic component of the deuteron wavefunction itself. Contrary 
to our previous analysis |13 the dd/dp ratio is no longer promising in this respect. The 
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next step in the clarification of this problem is clearly related to experimental data 
that e.g. can be taken at KEK. 



The authors like to thank A. A. Sibirtsev for valuable hints and inspiring discussions. 
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Figure Captions 



Fig. 1: The antiproton cross section as a function of the invariant energy above 
threshold ^/s — 4m. a) The open circles represent the results from the LUND string 
model (LSM) [15|] in comparison to the experimental data for the inclusive production 
[|i~4H (full squares). The crosses represent the results from the LSM for the exclusive 
channel pp —>■ pppp and the solid line is the result from the one-boson-exchange (OBE) 
model described in the text including p°, u, and tt° exchange while the dotted line 
stands for the 7r° exchange contribution, only. 

b) The phase-space model - including only the integrated 4-body phase space - is 
displayed by the dotted line and practically coincides with the result from the OBE 
model up to — 4m ^0.7 GeV. The solid line represents the fit from eq. (8) while 



the dashed line shows the approximation from Batko et al. |T7 



Fig. 2: The one-boson exchange model for the exclusive process pp — > pppp. 

Fig. 3: a) The Lorentz invariant differential cross section E d 3 a/dp 3 in (mb c 3 GeV~ 2 ) 
at 9i a b = 0° for the reaction d + p—^p + Xata deuteron momentum Pa = 10 GeV/c. 
The solid and dashed lines are the calculations with the Paris d.w.f. p4j using the 



parametrization of the elementary cross section from 17] and eq. (8), respectively. 



The dashed-dotted and dotted lines are the calculations with the Paris d.w.f. and a 



3% admixture of a non-nucleonic component using the parametrization from [I7| 
and eq. (8), respectively. 

b) Lorentz invariant differential cross section within the same units and notations as 
in a) for the reaction dd — * pX at P^= 7 GeV/c. 

Fig. 4: Ratio of Lorentz invariant antiproton spectra at 6i a b = 0° from d + d and 
d + p reactions at Pd =9.5 GeV/c for the Paris d.w.f. (dashed line) and the sum of the 
contribution from the Paris d.w.f. and a 3% admixture of a non-nucleonic component 



according to ref. [T3| . 
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